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Abstract. We discuss the tensor product construction for complexes of differential forms 
and show how it can be applied to define shape functions and degrees of freedom for finite 
element differential forms on cubes in n dimensions. These may be extended to curvilinear 
cubic elements, obtained as images of a reference cube under diffeomorphisms, by using the 
pullback transformation for differential forms to map the shape functions and degrees of 
freedom from the reference cube to the image finite element. This construction recovers 
and unifies several known finite element approximations in two and three dimensions. In 
this context, we study the approximation properties of the resulting finite element spaces 
in two particular cases: when the maps from the reference cube are affine, and when they 
are multilinear. In the former case the rate of convergence depends only on the degree of 
polynomials contained in the reference space of shape functions. In the latter case, the rate 
of approximation is degraded, with the loss more severe for differential forms of higher form 
degree. 



1. INTRODUCTION 

Finite element simulations are often performed using meshes of quadrilaterals in two di- 
mensions or of hexahedra in three. In simple geometries the meshes may consist entirely 
of squares, rectangles, or parallelograms, or their three-dimensional analogues, leading to 
simple data structures. However, for more general geometries, a larger class of quadrilateral 
or hexahedral elements is needed, a common choice being elements which are the images of 
a square or cube under an invert ible multilinear map. In such cases, the shape functions for 
the finite element space are defined on the square or cubic reference element and then trans- 
formed to the deformed physical element. In [5j it was shown in the case of scalar (H^) finite 
elements in two dimensions that, depending on the choice of reference shape function space, 
this procedure may result in a loss of accuracy in comparison with the accuracy achieved 
on meshes of squares. Thus, for example, the serendipity finite element space achieves only 
about one half the rate of approximation when applied on general quadrilateral meshes, as 
compared to what it achieves on meshes of squares. The results of [5j were extended to three 
dimensions in [15j. The case of vector (i/(div)) finite elements in two dimensions was studied 
in [6j . In that case the transformation of the shape functions must be done through the Piola 
transform and it turns out the same issue arises, but even more strongly, the requirement 
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on the reference shape functions needed to ensure optimal order approximation being more 
stringent. Some results were obtained for i7(curl) and i7(div) finite elements in 3-D in [14] . 

The setting of finite element exterior calculus (see [Hill [8]) provides a unified framework 
for the study of this problem. In this paper we discuss the construction of finite element 
subspaces of the domain HK^ of the exterior derivative acting on differential fc-forms, < 

< n, in any number n of dimensions. This includes the case of scalar finite elements 
(A; = 0), finite elements (in which the Jacobian determinant enters the transformation, 
k = n)^ and, in three-dimensions, finite elements in i/(curl) {k = 1) and H{div) (fc = 2). 

The paper begins with a brief review of the relevant concepts from differential forms on 
domains in Euclidean space. It then describes the tensor product construction for differential 
forms, and complexes of differential forms. These will be used to construct reference shape 
functions and degrees of freedom. First, in Section [4| we discuss the construction of finite 
element spaces of differential forms and the use of reference domains and mappings. Com- 
bining these constructions, in Section [5| we define the Q~A^ finite element spaces, which 
may be seen to be the most natural finite element subspace of HA^ using mapped cubical 
meshes. In Section |6] we obtain conditions for 0{h^^^) approximation in for spaces of 
differential fc-fo rms on mapped cubic meshes. More precisely, in the case of multilinear map- 
pings Theorem 6.1 shows that a sufficient condition for 0{h^^^) approximation in is that 
the reference finite element space contains Q~^^A^. Again we see loss of accuracy when the 
mappings are multilinear, in comparison to affine (for which VrA^ guarantees approximation 
of order r + 1), with the effect becoming more severe for larger k. In the final section we de- 
termine the extent of the loss of accuracy for the Q~A^ spaces. We discuss several examples 
of finite element spaces, including standard Qr spaces (fc = 0), VrA^ spaces, and serendipity 
spaces SrA^ which have been recently described in [3j. 



2. Preliminaries on differential forms 

We begin with a brief review of some basic notations, definitions, and properties of 
differential forms defined on a domain in M^. A differential fc-form on is a func- 
tion S Alt^M^, where Alt^M^ denotes the space of alternating fc-linear forms mapping 
(W)^ = W X • • • X W ^ R. By convention Alt^M"" = M, so diflFerential 0-forms are simply 
real- valued functions. The space Alt^ has dimension (^) for < /c < n and vanishes for 
other values oi k. A basis is formed by the basic alternators rfx^ := rfx^^ A • • • A rfx^'^, for 
a = (cTi, . . . , cT/e) belonging to the set S(A;, n) of increasing maps from {1, . . . , A;} to {1, . . . , n}. 
Here rfx' G Alt^ W = denotes the functional dx'^{v) = v"^ for v = {v^ , . . . , v"^) G and 

the wedge product of alternating forms is the skew part of the tensor product: 

/ A ^ = (^ ^ ^) skw(/ ®g), fe Alt' M", g G Alt' R\ 
Thus a general differential fc-form on S can be expressed uniquely as 
(1) f= J2 Udx''= J2 U-a„dx''' A-'-Adx'^'', 

for some coefficient functions on S. In particular, a differential n-form can be expressed 
diS f = gdx where the volume form dx = dx^ A • • • A dx^ and the coefficient is a function 
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on S. As long as the coefficient is integrable, the integral 

f = gdx 



of the differential form is defined and has the value the notation suggests. 

If J-'{S) is some space of real- valued functions on then we denote by TA^{S) the 
space of differential fc- forms with coefficients in J-'{S). This space is naturally isomorphic 
to J^{S) Alt^M^. Examples are the spaces C'^A^{S) of smooth fc- forms, L^A^(S') of 
fc-forms, H^A^{S) of /c- forms with coefficients in a Sobolev space, and VrA^{S) of forms with 
polynomial coefficients of degree at most r. The space L'^A^{S) is a Hilbert space with inner 
product 

and similarly for the space H^A^{S)^ using the usual Sobolev inner-product 

\a\<r 

where the sum is over multi-indices a of order at most r. 
The exterior derivative of the fc-form ([T]) is the {k + l)-form 

assuming that the indicated partial derivatives exist. We may view the exterior derivative 
as an unbounded operator L^A^(aS) L'^A^^^{S) with domain 

HA^{S) :={f e L^A^{S) \ df G L^A^^\S) }, 

which is a Hilbert space when equipped with the graph norm 

{fi9)HAk(s) = {f^9)L^Ak(S) + {df,dg)L2j^k+i(^s)' 

The complex 

(2) ^ HA'^iS) A HA\S) A . . . ^ HA'^iS) 

is the de Rham complex on S. 

A very important construction is the puUback of a differential form under a mapping. Let 
aS be a domain in and F = (F-^, . . . , F^) a mapping of S into a domain S in some 
M^. Given a differential fc-form v on aS, its puUback F^'v is a differential fc-form on S. If 

v= ^2 Vi^...i^dx'^ A • • • Adx'^, 

l<ii<"-<ifc<m 

then 

^ r)F^^ r)P^k 

(3) F*v= y y ivi,...i, o F)— —^dx^'^■■■^dx^K 

l<ii<---<ifc<m ji,...,ifc=l 
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The puUback operation satisfies {Go py = F* oG*, and, if F is a diffeomorphism, (F*)~^ = 
(F~^y. The puUback preserves the wedge product and the exterior derivative in the sense 
that 

F*(^ A ^) = F*^ A F*^, dF% = F^'dv. 

If F is a diffeomorphism of S onto S which is orientation preserving (i.e., its Jacobian 
determinant is positive), then the puUback preserves the integral as weU: 



Js Js 



An important apphcation of the puUback is to define the trace of a differential form on a 
lower dimensional subset. If aS is a domain in and / is a subset of the closure S and also 
an open subset of a hyperplane of dimension < n in M^, then the puUback of the inclusion 
map f '-^ S defines the trace operator tif taking fc- forms on S to /c-forms on /. 

Before continuing, we recall that vector proxies exist for differential 1-forms and (n — 1)- 
forms on a domain in M^. That is, we may identify the 1-form Vi dx^ with the vector field 
{vi^ . . . ^Vn) : S ^ M^, and similarly we may identify the (n — l)-form 

^^{-lyvi dx^ A"'dx' A dx"" 

i 

with the same vector field. A 0-form is a scalar function, and an n-form can be identified 
with its coefficient which is a scalar function. Under these identifications, the de Rham 
complex ([2]) on a three-dimensional domain becomes 

^ H\S) ^ //(curl, S) ^ //(div, S) ^ L\S) ^ 0. 

The puUback of a scalar function /, viewed as a 0-form, is just the composition: v{x) = 
'u(F(x)), X G S. If we identify scalar functions with n-forms, then the puUback is v{x) = 
det[D F {x)]v (^F (x)) where DF{x) is the Jacobian matrix of F at x. For a vector field 
viewed as a 1-form or an (n — l)-form, the puUback operation corresponds to 

v{x) = [DF{x)fv{F{x)), v{x) = adj[F>F(x)]^(F(x)), 

respectively. The adjugate matrix, adj A, is the transposed cofactor matrix, which is equal 
to {det A)A~^ in case A is invertible. The latter formula, representing the puUback of an 
(n — l)-form, is called the Piola transform. 

Next we consider how Sobolev norms transform under puUback. If F is a diffeomorphism 
of S onto S smooth up to the boundary, then each Sobolev norm of the puUback F^v can 
be bounded in terms of the corresponding Sobolev norm of v and bounds on the partial 
derivatives of F and F~^. Specifically, we have the following theorem. 

Theorem 2.1. Let r be a non-negative integer and M > 0. There exists a constant C 
depending only on r, M , and the dimension n, such that 

whenever S are domains in and F : S ^ S is a (7^+^ diffeomorphism satisfying 

max |F|^^(^)<M, |F-Vi(5) < M. 



l<s<r+l 
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Proof. From ([3]), 

\\F-v\\,.,.,,,= J: 5^ /|I."[K...,oF)— ...^|(«)|^<i4. 

l<ii<--<ik<m \a\<r ^ 

Using the Leibniz rule and the chain rule, we can bound the integrand by 

where C depends only on maxi<5<^+i \F\^s ^s) so can be bounded in terms of M. Chang- 
ing variables to x = Fx in the integral brings in a factor of the Jacobian determinant of 
which is also bounded in terms of M, and so gives the result. □ 

A simple case is when F is a dilation: F{x) = hx for some h > 0. Then ^ becomes 

{F%){x)= v^{hx)h^ dx"" , 

and therefore, 

D''{F%){x)= {D''v^){hx)h'^^ dx"" , 

where r = |q{|. We thus get the following theorem. 

Theorem 2.2. Let F be the dilation Fx = hx for some h> S a domain in , S = F{S), 
and a a multi-index of order r. Then 

3. Tensor products of complexes of differential forms 

In this section we discuss the tensor product operation on differential forms and complexes 
of differential forms. The tensor product of a differential fc-form on some domain and a 
differential /-form on a second domain may be naturally realized as a differential (/c + /)-form 
on the Cartesian product of the two domains. When this construction is combined with the 
standard construction of the tensor product of complexes, we are led to a realization of the 
tensor product of subcomplexes of the de Rham subcomplex on two domains as a subcomplex 
of the de Rham complex on the Cartesian product of the domains (see also P[ [10]). 

We begin by identifying the tensor product of algebraic forms on Euclidean spaces. For 
m, n > 1, consider the Euclidean space ]R^+^ with coordinates denoted (xi, . . . , x^, yi, . . . , yn) 
The projection tti : ]R^+^ on the first m coordinates defines, by puUback, an injec- 

tion nl : Alt^M^ Alt^]R^+^, with the embedding of tt* : Alt^M^ Alt^]R^+^ defined 
similarly. Therefore we may define a bilinear map 

Alt^ X Alt^ ^ Alt^+^ M"^+", (/X, u) ^ 7r*/x A tt^z/, 

or, equivalent ly, a linear map 

Alt^ W Alt^ W Alt^+^ ]R^+^, fi^u^nlfiA n^u. 

This map is an injection. Indeed, a basis for Alt^M^ (8) Alt^M^ consists of the tensors 
dx^ dy^ with a G m) and r G n), which simply maps to dx^ A dy^ ^ an element of 
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the standard basis of Alt^+^ M^+^. In view of this injection, we may view the tensor product 
of /i z/ of a /c-form /x on and an /-form on as a (A; + /)-form on ]R^+^ (namely, we 
identify it with tt^/x A 7r2Z^). 

Next we turn to differential forms defined on domains in Euclidean space. Let u he Sl 
differential fc-form on a domain S C and v a differential /-form on T C M^. We may 
identify the tensor product u(^v with the differential {k + /)-form n^u A n^v on S xT where 
TTs S X T ^ S and ttt ' S x T ^ T are the canonical projections. In coordinates, this 
identification is 

{J2 U dx^) ® 9r dy^ = Y.f-®3r dx^ A dy\ 



Note that the exterior derivative of the tensor product is 

dsxriTTsU A Ti^v) = dsxriTTsu) A n^v + (-l)^7rj2i A rf^xT(vr^^) 

(4) , 

= 7Ts{dsu) A TT^V + {-IfTT^^U A 7T^{dTv). 

Having defined the tensor product of differential forms, we next turn to the tensor product 
of complexes of differential forms. A subcomplex of the de Rham complex ([2]), 

(5) ^ y° 4 4 • • • 4 ^ 0, 

is called a de Rham subcomplex on S. This means that, for each fc, C HA^{S) and d 
maps into 1/^+^. Suppose we are given such a de Rham subcomplex of S and also a 
de Rham subcomplex, 

(6) O^W^ ^W^ ^ Aw ^0, 

on T. The tensor product of the complexes ^ and ^ is the complex 

(7) ^ {V ®wf A {V ®wY A • • • A {V ® ^ o, 

where 

(8) {V W)^ := ^{V'^W), A; = 0, . . . , m + n. 

i-\-j=k 

and the differential (V W)^ {V (S) W)^^'^ is defined by 

d{u 0v) = d'gu ®v + {-lyu d^r^v, u ^V\v ^ 

In view of the identification of the tensor product of differential forms with differential forms 
on the Cartesian product, the space {V ®WY in ^ consists of differential fc-forms on x T, 
and, in view of Q, the differential is the restriction of the exterior derivative on HA^{S x T). 
Thus the tensor product complex ([7| is a subcomplex of the de Rham complex on x T. 

4. Finite elements 

Let be a bounded domain in M^, Q < k < n. As in [llj, a finite element space of fc-forms 
on Vt is assembled from several ingredients: a triangulation T of f^, whose elements we call 
finite elements^ and, for each finite element a space V{K) of shape functions on and 
a set of degrees of freedom. We now describe these ingredients more precisely. 

For the triangulation, we allow the finite elements to be curvilinear polytopes. This means 
that each K is the image Fk{K) of an n-dimension polytope K (so a closed polygon in two 
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dimensions and a closed polyhedron in three dimensions) under a smooth invertible map Fk 
of K into R^. The faces of K are defined as the images of the faces of and the requirement 
that T be a triangulation of Vt means that Vt = UxgT ^ ^^^^ intersection of any 
two elements of T is either empty or is a common face of both of some dimension. 

The space V{K) of shape functions is a finite-dimensional space of /c- forms on K. The 
degrees of freedom are a unisolvent set of functionals on V{K), or, otherwise put, '^{K) is a 
basis for the dual space V{Ky. Further, each degree of freedom is associated to a specific 
face of X, and when two distinct elements Ki and K2 intersect in a common face /, the 
degrees of freedom of Ki and K2 associated to the face / are in 1-to-l correspondence. 

The finite element functions are then defined as the differential forms on Q which belong 
to the shape functions spaces piecewise, and for which corresponding degrees of freedom 
are single- valued. That is, to define a finite element function u we specify, for all elements 
shape functions Uk ^ ^(^) satisfying ^i{uki) = ^2(^X2) whenever Ki and K2 share 
a common face and G ^i{K) and ^2 ^ ^2{K) are corresponding degrees of freedom 
associated to the face. Then u is defined almost everywhere by setting its restriction to the 
interior of each element K to be uk- The finite element space V{T) is defined to be the 
space of all such finite element functions. 

The degrees of freedom associated to faces of dimension < n determine the inter-element 
continuity imposed on the finite element functions. Specifying the continuity in this way 
leads to finite element spaces which can be efficiently implemented. We note that the finite 
element space is unchanged if we use a different set of degrees of freedom, as long as the 
span of the degrees of freedom associated to each face is unchanged. Thus we shall usually 
specify these spans, rather than a particular choice of basis for them. 

When a finite element K is presented as Fk{K) using a reference element K dW^ and 
a diffeomorphic mapping Fx . K ^ K (as for curvilinear polytopes), it is often convenient 
to specify the shape functions and degrees of freedom in terms of the reference element and 
the mapping. That is, we specify a finite dimensional space V{K) of differential A;-forms on 
the reference shape functions^ and define V{K) — {F^^yV{K)^ the puUback under the 
diffeomorphism F^^ \ K ^ K. Similarly, given a space 5(^) of degrees of freedom for V{K) 
on K and a reference degree of freedom ^ G associated to a face / of we define ^ 

by ^{v) = C{F^v) for V G V{K) and associate ^ to the face / = F^if) of K. In this way we 
determine the degrees of freedom In the present paper we shall be concerned with the 

case where each element of the mesh is presented as the image of the unit cube K = [0^1]^ 
under a multilinear map (i.e., each component of Fk is a polynomial of degree at most one 
in each of the n variables). A special case is when Fk is affine, in which case the elements 
K are arbitrary parallelotopes (the generalization of parallelograms in two dimensions and 
parallelepipeds in three dimensions). In the general case, the elements are arbitrary convex 
quadrilaterals in two dimensions, but in more dimensions they may be truly curvilinear: the 
edges are straight, but the faces of dimension two need not be planar. 

We need some measure of the shape regularity of an element K. Since the shape regularity 
should depend on the shape, but not the size, of the element, we require it to be invariant 
under dilation. Therefore, let Hk = diam(i^^) and set K = h^K^ which is of unit diameter. 
Assuming only that the element is a Lipschitz domain, the Bramble-Hilbert lemma ensures 
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that for r > integer there exists a constant C such that 

(9) inf _ 11^ - p\\l2^k) < u G H'^^\K). 

peVr(K) 

We define C{K^r) to be the least such constant, and take this as our shape regularity 
measure. If the domain K is convex then C{K^r) can be bounded in terms only of r and 
the dimension n; see [18j. If the domain is star-shaped with respect to a ball of diameter at 
least Shx for some 5 > 0, then C{K^ r) can be bounded in terms of r, n, and 5] see [13]. 

The next theorem will play a fundamental role in proving approximation properties of 
finite element differential forms. 

Theorem 4.1. Let K he a hounded Lipschitz domain in with diameter hx and let r he a 
non-negative integer. Then there exists a constant C depending only on r, and the shape 
constant C{K^r), such that 

inf \\u - p\\L'^Ak(K\ < Ch'j^'^\u\Hr+i/^k(x), u G H^'^^K^iK), 

for < k < n. 

Proof. Let K = h'^K, lei F \ K ^ K he the dilation, and set u = F*2i. Applying ^ to 
each component of we obtain p G Vr^^{K) with 

(10) ||F*2i - pWl^K^^k) < ClF^'ulHr+iAk^Ky 

Set q = {F~'^yp. Since F~^ is a dilation, q G VrA^{K), and clearly F"'{u - q) = F^'u - p. 



Combining (|10[) and Theorem 2.2 gives the desired estimate. □ 



5. Tensor product finite element differential forms on cubes 

In this section, we describe the construction of specific finite element spaces of diflFerential 
forms. We suppose that the triangulation consists of curvilinear cubes, so that each element 
K of the triangulation T is specified by a smooth diflFeomorphism F^ taking the unit cube 
K = [0, 1]^ onto K C M^. As described above, the shape functions and degrees of freedom 
on K can then be determined by specifying a space V{K) of reference shape functions and 
a set S(i^) of reference degrees of freedom. Now we apply the tensor product construction 
of Section [3] to construct the reference shape functions and degrees of freedom. 

Fix an integer r > 1. With / = [0, 1] the unit interval, let P^A°(/) = Vr{I) denote 
the space of polynomial functions (0- forms) on / of degree at most r, and let Vr-iA.^{I) = 
Vr-i{I) dx denote the space of polynomial 1-forms of degree at most r — 1. Connecting these 
two spaces with the exterior derivative, we obtain a subcomplex of the de Rham complex on 
the unit interval: 

(11) Vr{I)^Vr-l{I)dx. 

Taking the tensor product of this complex with itself, we obtain a subcomplex of the de Rham 
complex on whose spaces we denote by Q~A^ = Q~A^{P): 

(12) Q;A' A Q;A' A Q;A^. 

Specifically, writing Vr^s = ^r,s(^^) for ^r(^) ® ^s{I)) we have on 

Q;A^ = Vr^r^ Q;A^ = Vr-i,r dx Vr^r-i dy , A^ = Vr-i,r-i dx A dy. 
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The first space Q~ is the tensor product polynomial space traditionally denoted Q^, and 
the last space is Qr-i dx A dy. 

To get spaces on the 3-D cube, we may further take the tensor product of the 2-D complex 



(12) with the 1-D complex (11), or, equivalently, take the tensor product of three copies of 



the 1-D complex, to obtain a de Rham subcomplex on the unit cube P: 
Here Q;A^ = Qr, Q'A^ = Q^-i dxAdyA dz, and 

Q~A^ = Vr-l,r,r dx 7^r,r-l,r dy © 7^r,r,r-l dz, 

Q~A^ = Vr,r-i,r-i dy A dz Vr-i,r,r-i dx A dz © Vr-i^r-i,r dx A dy. 

The extension of this construction to higher dimensions is clear, yielding spaces Q~A^{I^) 
for < < n, n > 1, spanned by quantities pdx^^ where a G n) and p is a polynomial 
of degree at most r in all variables and of degree at most r — 1 in the variables x^\ More 
precisely. 



Q;A\r)= 



1=1 



dx""' A'" Adx^'K 



where 



1, i e {ai, . . . ,ak}, 
0, otherwise. 



In the case = 0, this space is understood to be 

n 
i=l 

i.e., the space conventionally referred to as Q^(/^). This definition makes sense also if r = 0, 
so QqA^Ii"^) = Qoil"^) = M is the space of constant functions. For A; > 0, QoA^{I'^) = 0. 
We also interpret this in the case n = 0, so is a single point. Then Q~A^(/°) is understood 
to be the space M of constants. It is easy to check that 

dim Q;A\r) = Q (r + l)"-^r^ 

in all cases < A; < n, r > 0. 

Let us now characterize the degrees of freedom of Q~A^{I^). The space VrA^{I) = Vr{I) 
has dimension r + 1. It has one degree of freedom associated to each of the two vertices 
p = and p = 1, namely the evaluation functional v ^ v{p). The remaining r — 1 degrees 
of freedom are associated to / itself, and are given by 

V J v{x)q{x) dx, q G Vr-2{I)- 

The space Vr-iA^{I) = Vr-i{I) has dimension r and its degrees of freedom are given by 

V v{x)q{x) dx, q G Vr-i{I)- 



I 
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The tensor product construction then yields degrees of freedom for Q~A^{I^). Namely, each 
degree of freedom for Q~A^{I^) is of the form Ci ® • • • ® Cn, and is associated to the Cartesian 
product face /i x . . . x where G is associated to the face fj of the interval I W j. 
Hence, a set of unisolvent degrees of freedom for Q~A^{I^) (r > 1, < A; < n) is given by 

(13) [ tTfv{x)Aq{x), q G Q;_,A'-'{f)^ 

for each face / of of degree d > k. 

Let us count that we have supplied the correct number of degrees of freedom. Since the 
number of faces of dimension d of is 2^~^(^), the total number of degrees of freedom is 

Substituting (^) (^) = (^) (^I^) and then changing the summation index from d to m = n — d 
we get: 

t (:) (: : 3 - - (:) e c ; ') ^-c- - d"-'- 

d=k \ / \ / \ / ^ / 

;^r^(r - 1 + 2)^-^ = dim Q;A^{r), 

as desired. 

Figure [l| taken from [2], shows degree of freedom diagrams for the Q~A^ spaces in two 
and three dimensions. In such diagrams, the number of symbols drawn in the interior of a 
face is equal to the number of degrees of freedom associated to the face. 

Now suppose that a mesh T of curvilinear cubes K = Fk{K) is given. We assume that 
whenever two elements Ki and K2 meet in a common face, say / = = i^Xsl/s)? 

then the map F^^ o i^xjj^ mapping /i onto /2 is linear. The assembled finite element space 
Q~ A^{T) consists of functions u whose restriction uk to K belongs to [F^^Y Q~ A^ [I^) and 
for which the corresponding degrees of freedom are single- valued on faces. The choice of 
degrees of freedom implies that if Ki and K2 share a face / of dimension > fc, then the 
traces ti fUx^ and tr/t^^s coincide. This is exactly the condition needed for u to belong 
to the space HA^iQ)^ i.e., for the exterior derivative of u to belong to L'^A^^^iyt) (see jH 
Lemma 5.1]). From the commutativity of the puUback with the exterior derivative, we see 
that dQ~ A^{T) C Q~A^+-^(T), and so we obtain a subcomplex of the de Rham complex 
on Vt\ 

Q;A\T) 4 Q;A\T) 4 ... 4 q;a^{T). 

The assembled finite element spaces Q~A^{T) are well-known, especially when the maps 
Fk are composed of dilations and translations, so the mesh consists of cubes. The space 
Q~A^{T) is the usual Qr approximation of H^{Q) and the space Q~A^{T) is the discon- 
tinuous Qr-i approximation of L'^{Q). In three dimensions, if we identify HA^ and HA^ 
with i7(curl) and i7(div), respectively, then Q~A^{T) and Q~A^{T) are the Nedelec edge 
element and face element spaces of the first kind, respectively [17] . 





Figure 1. Degree of freedom diagrams for spaces in two and three dimensions. 
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6. Approximation properties on curvilinear cubes 

We now return to the situation of Section [4] and consider the approximation properties 
afforded by finite element spaces of diflFerential /c- forms on meshes of curvihnear cubes. Thus 
we suppose given: 

(1) A space of reference shape fc-forms V{K) and a set of reference degrees of freedom 

on the unit cube K; 

(2) A mesh T of the domain 

(3) For each X G T a diffeomorphism of K onto K. 

As described in Section [4| these determine a space of fc-form shape functions V{K) on each 
element K and an assembled finite element space V{T) consisting of /c- forms on Q which 
belong piecewise to the V{K ). If it happens that V{K) contains the full polynomial space 



VrA^{K), then Theorem 4.1 gives us the approximation result 



(14) inf \\u - w||L2Afc(x) < Ch'j^^lulHr+iAk^K), 

with a constant C depending only on r, n, and the shape constant C{K^r). 



We now show how to extend estimate (14) to the global space V{T). More precisely, under 
the same condition that for all X G T the space V{K) contains the full polynomial space 
VrA.^{K)^ we want to show that 

(15) inf ll'U — '^ILsAfc ^ Ch'^^^\u\Hr+iAk^ 

veviT) 

where h is as usual the maximum of the /i^'s and the constant C depends only on r, n, and 
the shape constants C{K^r). This result can be obtained by extending to this setting the 
classical construction of the Clement interpolant presented in [12]. An analogous extension 
has been presented in [7, Sec. 5.4] in the case of simplicial meshes. A Clement-like operator 



Tih ' L^A^ ^{^) can be defined as in (13), where tifV^x) is replaced by the trace of the 



projection of v onto Vr^^{S)^ S being the union of elements in T containing the face 



/. Following [12], crucial conditions in order to get (15) are that global polynomial forms 
are preserved by tt/^, that is TThW = w whenever w belongs to P^A^(f]), and suitable scaling 
estimates which in our case are consequences of Theorems |2.1| and |2.2[ 

Thus we are led to ask what conditions on the reference shape functions V{K) and the 
mappings Fx ensure that V{K) = {F-^yV{K) contains VrA^{K). The following result, 
which is straightforward to prove at this point, answers this question for afl&ne and multilinear 
maps. 

Theorem 6.1. Suppose that either 

(1) Fk is an affine diffeomorphism and that VrA.^{K) C V{K); or 

(2) Fk is a multilinear diffeomorphism and Q~^j^A^{K) C V{K) 



Then V{K) contains VrA^{K) and so ^ and ^ hold. 



Proof Let us write F for Fx- The requirement that V{K) = (F"^)*1/(K) contains P^A^(X) 
is equivalent to requiring that F*(7^^A^(K)) C V{k). If F is affine, then F*(7^^A^(K)) C 
VrA.^{K)^ as is clear from ([3]). The sufficiency of the first condition follows. 



TENSOR PRODUCT FINITE ELEMENT DIFFERENTIAL FORMS 



13 



For F multilinear, we wish to show that F*(P^A^(X)) C A^(i^), for which it suffices 

to show that F*{pdx'') G Q;^j,A^{K) if p G Vr{K) and a G T.{k,n). From g it suffices to 
show that 

(16) " • • • ^^'^ ^ • • • ^ ^^'^ ^ 2;+.A'(^). 

Now p G and F multilinear imply that poF G Qr{K). Moreover, dF'^^ /dx^^ is multilinear, 
but also independent of x^K Therefore the product, [p o F)dF'^^ / dx^^ • • - dF'^^ / dx^^ ^ is of 
degree at most r + in all variables, but of degree at most r + — 1 in the variables x^K 
Referring to the description of the spaces Q~ spaces derived in Section [5| we verify (16). □ 



Note that the requirement on the reference space V{K) to obtain 0{h^^^) convergence 
is much more stringent when the maps Fk are only assumed to be multilinear, than in the 
case when they are restricted to be affine. Instead of just requiring that V{K) contain the 
polynomial space VrA^{K)^ it must contain the much larger space Q~^^A^(K). For 0-forms, 

the requirement reduces to inclusion of the space Qr{K). This result was obtained previously 
in [5] in two dimensions, and in p^, in three dimensions. The requirement becomes even 
more stringent as the form degree, fc, is increased. 

7. Application to specific finite element spaces 



Theorem 6.1 gives conditions on the space V[K) of reference shape functions which ensure 
a desired rate of approximation accuracy for the assembled finite element space V{T). In 
this section we consider several choices for V{K) and determine the resulting implied rates 
of approximation. According to the theorem, the result will be different for parallelotope 
meshes, in which each element is an affine image of the cube X, and for the more general 
situation of curvilinear cubic meshes in which element is a multilinear image of K (which we 
shall simply refer to as curvilinear meshes in the remainder the section). Specifically, if s is 
the largest integer such that 

(17) Vs-iiy^k) c v{k), 

then the theorem implies a rate of s (that is, error bounded by 0{h^)) on parallelotope 
meshes, while if s is the largest integer for which the more stringent condition 

(18) Q-^u-i^\k) c v{k), 

holds, then the theorem implies a rate of s more generally on curvilinear meshes. 

7.1. The space Q^. First we consider the case of 0-forms {H^ finite elements) with V{K) = 
Q~hP{K)^ which is the usual Qr{K) space, consisting of polynomials of degree at most r 



in each variable. Then (17) holds for 5 = r + 1, but not larger. Therefore we find the 



approximation rate to be r + 1 on parallelotope meshes, as is well-known. Since = 0, 



(18) also holds for 5 = r + 1, and so on curvilinear meshes we obtain the same rate r + 1 
of approximation. Thus, in this case, the generalization from parallelotope to curvilinear 
meshes entails no loss of accuracy. 
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7.2. The space A^. The case of n-forms with V{K) = Q^A^{K) is quite different. If 



we identify n-forms with scalar-valued functions, this is the space Qr-i{K). Hence (17) 



holds with 5 = r, and so we achieve approximation order 0{h^) on parallelotope meshes. 



However, for A; = n, (18) holds only if5<r — n + 1, and so Theorem 6.1 only gives a rate 
of r — n + 1 on curvilinear meshes in this case, and no convergence at all if r < n — 1. Thus 
curvilinear meshes entail a loss of one order of accuracy compared to parallelotope meshes 
in two dimensions, a loss of two orders in three dimensions, etc. This may seem surprising, 
since if we identify an n-form v dx^ A • • • A dx^ with the function which is a 0-form, the 
space Q~A^{K) corresponds to the space Qr-i{K) = Q~_iA^{K)^ which achieves the same 
rate r on both classes of meshes. The reason is that the transformation of an n-form on K 
to an n-form on K is very different than the transform of a 0-form. The latter is simply 
V = V o F while the former is 

vdx^ A'" Ax"" = {vo F)(det F) dx^ A " - A dx"". 
7.3. The spaces Q~ A^. Next we consider the case V{K) = Q~ A^{K) for k strictly between 



and n. Then (17) holds with s = r and (18) holds with s = r — k -\- 1 (but no larger). 



Thus, the rates of approximation achieved are r on parallelotope meshes and r — k -\- 1 for 
curvilinear meshes. In particular, if r < fc, then the space provides no approximation in 
the curvilinear case. For example, in three dimensions the space Q^A^(X), which, in con- 
ventional finite element language is the trilinearly mapped Raviart-Thomas-Nedelec space, 
affords no approximation in L^. This fact was noted and discussed in [16]. 

7.4. The space VrA^. Since the space HA^{Vt) does not require inter-element continuity, 
the shape functions V{K) = VrA^{K) may be chosen as reference shape functions (with all 



degrees of freedom in the interior of the element). In this case, (17) clearly holds for 5 = r + 1, 
giving the expected 0{h^^^) convergence on parallelotope meshes. But (18) holds if and only 
if QsJ^n-2{K) C Vr{K)^ which holds if and only if n{s + n — 2) < r. In two dimensions this 
condition becomes 5 < r/2, and so we only obtain the rate of [r/2j , and no approximation at 
all for r = or 1. In three dimensions, the corresponding rate is 0[r/3j — 1, requiring r > 6 
for first order convergence on curvilinear meshes. For general n, the rate is [r/nj — n + 2. 

7.5. The serendipity space 5^. The serendipity space 5^, r > 1, is a finite element sub- 
space of (i.e., a space of 0- forms). In two dimensions and, for small r, in three dimensions, 
the space has been used for many decades. In 2011, a uniform definition was given for all 
dimensions n and all degrees r > 1 [3j. According to this definition, the shape function 
space Sr{K) consists of all polynomials of superlinear degree less than or equal to r, i.e., for 
which each monomial has degree at most r ignoring those variables which enter the mono- 
mial linearly (example: the monomial x^yz^ has superlinear degree 5). From this definition, 
it is easy to see that Vr{K) C Sr{K) but Vr+i{K) ^ Sr{k). Thus, from (|l7|), the rate of 
approximation of Sr on parallelotope meshes is r + 1. Now Qi{K) C Sr{K) for all r > 1, 
but for 5 > 2, Qs{K) contains the element (xi • • - XnY of superlinear degree ns. It follows 



that (18) holds if and only ii s < max(2, [r/nJ + 1), which gives the L rate of convergence 



of the serendipity elements on curvilinear meshes. This was shown in two dimensions in 
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7.6. The space cS^A^. In [4J, Arnold and Awanou defined shape functions and degrees 
of freedom for a finite element space 5^A^ on cubic meshes in n-dimensions for all form 
degrees k between and n, and polynomial degrees r > 1. The shape function space 
SrA.^{K) they defined contains P^A^(X), so the assembled finite element space aflFords a 
rate of approximation r + 1 on parallelotope meshes. In the case = n, SrA.^{K) in fact 
coincides with P^A^(i^), so, as discussed above, the rate is reduced to [r/nj — n + 2 on 
curvilinear meshes. In the case = 0, 5^A° coincides with the serendipity space 5^, and 
so the rate is reduced to max(2, [r/nJ + 1) on curvilinear meshes in that case. In n = 2 
dimensions, this leaves the space 5^A^, which is the BDM space on squares, for which the 
reference shape functions are Vr^^{K) together with the span of the two forms d{x^^^y) 

) we note that Q~K^[K) contains x^~^y^dx^ which 
is contained in SrA^{K) only if 2^ — 1 < r. Thus the rate of approximation of 5^A^ on 
curvilinear meshes in two dimension is [(r + l)/2j. 
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